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A Critical Review of Charged Particle Astronomy at Saturn: 
The Evidence for Co-orbiting Material in the Inner Satellite System 


SUMMARY 


,Y.«^ 

IN - 

Agreement NAG 2-528 was designed to support Dr. John F. Cooper's work on / 
reassessment of die charged particle observations from Pioneer and Voyager at Saturn with 
a view towards providing limits on the amount of "unseen" dust and debris that may exist in 
the Saturnian system. Such estimates are crucial for planning the Cassini tour of Saturn. 

The data from Pioneer- 1 1 and Voyager were reviewed, intercompared and correlated with 
model predictions to set limits on the matter distribution. The results were presented at 
several meetings and incorporated into the Cassini Proposal Information Package. 
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TECHNICAL DISCUSSION 


prepared by: John F. Cooper 

Department of Physics and Astronomy 
Louisiana State University 
Baton Rouge, Louisiana 70803 

1. Introduction 

The planetary system and magnetosphere of Saturn has been studied in-situ by three 
flyby missions: Pioneer 11 in 1979 [Opp, 1980], Voyager 1 in 1980 [Stone and Miner, 1981], 
and Voyager 2 in 1981 [Stone and Miner, 1982]. All three spacecraft were fully instrumented 
for measurements of energetic charged particles, but only Pioneer 11 covered the full range 
of the magnetosphere into the inner planetary system as shown in Figure 1. The first 
measurements of high energy trapped radiation were reported by the Pioneer 11 investigators 
[Trillins et al., 1 980 j F dim s and Mcllwain, 1980 , Simpson et o/., 1980a, b, 1981, Van Allen et 
al., 1980a,b] and further measurements with higher resolution in energy and direction were 
made by Voyager 2 [Krimigis et al., 1982, 1983; Krimigis and Armstrong, 1982; McDonald 
et al., 1980; Schardt and McDonald, 1983; Trainor et al., 1980; Vogt et al., 1982] The nearly 
equatorial flyby trajectory of Pioneer 11 was especially well-suited for measurements of 
radial distributions and time dependence (i.e., inbound versus outbound) of stably trapped 
particles, while the Voyager 2 trajectory covered a larger range of magnetic latitudes at 
L > 2.73 and was more useful for measurements of pitch angle distributions [e.g., Schardt 
and McDonald, 1983; Schardt, 1983]. . - r 

Ongoing plans [Stetson, 1989; Nicholson, 1989] for the next mission to Saturn, Cassini, 
call for an four-year tour of the Saturn system by the Cassini orbiter spacecraft after re- 
lease of the Titan probe. If properly instrumented, the orbiter tour will allow very detailed 
surveys of trapped particles within the magnetosphere and over the main rings. While pro- 
viding valuable new scientific information on the magnetospheric trapped radiation and its 
interaction with rings, satellites, other co-orbiting bodies, and dust in the Saturn system, 
charged particle measurements may also provide a critical capability for assessing proba- 
ble abundances, or upper limits thereof, for material posing potentially significant impact 
hazards. Since the tour will require many crossings of the equatorial ring plane within and 
inwards of the E-ring region, there is a clear need to evaluate the relative merits of different 
radial positions for these crossings, so that that impact hazards may be minimized. 

The present review specifically addresses the most relevant observations of energetic 
trapped particles from Pioneer and Voyager for the indirect survey of co-orbiting material 
via trapped particle absorption. While absorption effects of the major rings and satellites 
are clearly evident in the observations, the most critical questions for Cassini concern limits 
that may be placed on abundances of diffuse dust and larger debris which are not well- 
constrained by other observations. Although Voyager and/or earth-based observations have 
measured the total optical opacities of the G- and E-rings, these opacities are dominated 
by contributions from small, micron-sized dust particles when measured with backscattered 
sunlight. Particles above a few hundred microns in size pose the greatest impact hazards, 
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make little contribution to optical opacity, and may be better assessed by absorption of 

trapped radiation. _ 

The most useful trapped particle observations for the present review fall in three cate- 
gories: (1) high energy protons, which have very long trapping times (i.e., years) and cor- 
respondingly high probabilities for interaction with low-opacity dust material (if present), 
(2) high intensity electrons at energies near a few MeV, which show significant response 
to magnetospheric perturbations on hour time scales and can only be used for sampling of 
high-opacity absorbers, and (3) higher energy electrons which may arise in part as secondary 
products of proton-dust interactions. The principal focii of discussion are (1), which has a 
reasonably well-known source rate from cosmic ray albedo neutron decay (CRAND), and 
the potential evidence for the existence of (3), which might provide the best constraints on 
dust abundances if measured by appropriate instruments on the Cassini orbiter. We will also 
discuss ” microsignature” observations which may indicate the presence of large abundances 
of localized (but unseen) dust clouds co-orbiting with the principal satellites, although uncer- 
tainties (viz-a-viz electron background) remain in interpretation of the charged particle data 
for some microsignatures. All of the above categories are recommended for measurements 
by orbiter instrumentation during the early phases of the Cassini mission. 

Despite the indirect nature of the charged particle measurements for detection of absorb- 
ing dust and/or larger bodies, such measurements offer an important advantage over direct, 
in-situ measurements with dust detectors. The charged particle surveys can be carried out 
rapidly during the first few passes of the orbiter through the inner magnetosphere over a 
wide radial range in relative safety away from the equatorial region. In contrast the in-situ 
measurements can only be made at selected equatorial crossings where one has already es- 
timated that the dust densities will be minimal. A complementary arrangement would be 
that the charged particle data be used for the initial wide-ranging surveys, while later in-situ 
observations would be planned for ring- plane crossings at specific points of interest. Regions 
of hi gh trapped proton intensities can be chosen for the initial aim points at the ring plane 
before and after SOI to minimize risks during the first surveys. On the other hand there may 
be a an understandable desire by mission planners to avoid any regions of interest for direct 
dust measurements that appear even marginally hazardous from the charged particle survey. 
Thus the charged particle data may have the more crucial role for hazard assessment. 

High resolution science data on species, energy spectra, and directionality of high energy 
trapped radiation would allow prolonged study of sources, radial transport, and sinks while 
also providing critical remote-sensing information for assessment of spacecraft hazards asso- 
ciated with selection of sites for ring plane crossings during the four-year mission. Although 
high energy measurements of trapped radiation have been made for many years in the ter- 
restrial magnetosphere, Saturn offers a unique new environment to study interactions of this 
radiation with the planetary environment, such that significant new contributions may be 
made to magnetospheric and planetary science. By the time the Cassin orbiter finishes its 
presently conceived mission at Saturn by about the year 2010, the earth s man-made rings 
wil be dense enough to significantly absorb CRAND protons in the earth’s magnetosphere 
[Konradi, 1988], so that the study of dust-magnetosphere interactions may become increas- 
ingly important for assessing hazards from debris impacts with respect to unmanned and 
manned space operations near Barth as well as at Saturn. 
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Figure 1: Pioneer 11 trajectory in the inner Saturn system [Simpson ei oi., 1 ] 



in c inner masnctospherc from Pioneer 11 for electrons 

Figure 2: Inbound cou “ t “ g ^ “ J 1980 ]. Directional components are shown for electrons 

> 0.45 MeV and protons > 80 MeV jFdh^ ^ ^ fidd action. The proton rates are averaged 

Tk< MW UW«S ” 1879-S2" U . »*">■*“*»“ to U. o.totol «P°» • 


Janus and Epimetheus. 
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2. Observations of Energetic Trapped Radiation 


Sources 

High energy trapped protons and electrons in the inner magnetosphere of Saturn have 
complex interactions with planetary rings, satellites, and dust in the Saturn System. Nuclear 
interactions of very high energy (> 20 GeV/N) cosmic ray protons and nuclei with Saturn’s 
main rings produce secondary particles at sub-GeV energies with high efficiency [Blake et al., 
1983; Cooper, 1983; Cooper et al, 1985]. The charged secondaries have low level intensities 
within the ring region due to rapid reabsorption by the ring material, while the neutral 
secondary components (neutrons and gamma rays) can radiate away from the rings. Beta 
decay of the secondary neurons (CRAND) provides the only known source of protons above 
10 MeV in the trapped radiation belts of the inner magnetosphere. 

The beta-decay protons are trapped with essentially all of the parent neutron energy and 
ran undergo further interactions with satellites and more diffuse material orbiting around 
Saturn near the equatorial ring plane. A significant fraction of the trapped protons are suf- 
ficiently energetic to produce pions from interi- ;tior. nth dust grains, thereby populating 
the radiation belts with highly relativistic (10-100 McV) electrons from pion-muon decays. 
Inward radial diffusion and betatron acceleration of sub-MeV electrons from the outer mag- 
netosphere provides a much higher intensity population of electrons at energies below five 
MeV. The relative intensities of these trapped secondary interaction products are in turn 
controlled by lifetimes for radial transport within the inner magnetosphere and re-absorption 
by satellites and diffuse material. 

Radial Profiles 

The Pioneer 11 data of Fillius et al. (1980) are shown in Figures 2 and 3 to illustrate 
the entirely different radial profiles of the trapped electron and CRAND proton components 
in the inner magnetosphere. Despite the very low injection rates from CRAND, the fluxes of 
Fjgli energy protons build up to high levels in radial zones not swept by the mam rings, the 
F-ring, the co-orbiting satellites Epimetheus and Janus, and the innermost principal moons, 
Mimas and Enceladus. Since the rate of source injection varies slowly across the proton 
trapping regions, the relative magnitudes of proton intensities in these regions are primarily 
due to diffusive effects [see Section 3]. There is no obvious effect of absorption by diffuse 
dust other than the inflection in the proton profile in the vicinity of the G ring near 2.8 Rs, 
although dust losses might be masked by radial variations arising from diffusion. 

In sharp contrast, the MeV electron profile in Figure 2 shows no long-term effect of 
satellite absorption and is affected only at the inner edge of the trapped radiation region 
in the F-ring region before complete absorption occurs in the main rings. The 1979S2 
microsignature [Fillius et al, 1980; Simpson et al., ixfi 0a,b; Van Allen et al, 1980a, b] near 
2.5 planetary radii (1 R s = 60,000 km) was a highly localized signature of absorption 
resulting from a probable near encounter with one of the co-orbiting satellites [e.g., Van Allen, 
1984]. The absorption in the F-ring region is most likely due to multiple close encounters 
with distributed absorbers in that region [e.g., Cuzzi and Burns, 1988]. Thus the observed 
electron profiles are affected primarily by localized, strong absorbers. 

Inward diffusion and acceleration produces a ’’inner core” population of MeV electrons 
in the inner magnetosphere which has distinctly different radial profiles as compared to 
the seed population in the outer magnetosphere at lower energies. Figure 4 shows radial 
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profiles for integral counting rates of electrons above 0.040 and 0.56 MeV from the P.oneer 
11 encounter (Van Allen et « 1, 1980a; Van Allen, 1984]. While the e ectrons at 0.040-0.56 
MeV show time dependent behavior outside the orbit of Tethys, the higher energy electrons 
show a monotonic increase inwards which is characteristic of sternly ,nw„d acceleration for 
electrons counted by detectors with fixed energy thresholds. Just outside Tethys s or i , 
two profiles converge (note the different vertical scales) into a sing e core P ro 
inwards which is dominated by electrons above the higher energy threshold. The more i rapid 
increases in the counting rates inside the orbit of Enceladus are unexplained but migh be due 
to changes in the electron spectra, due to preferential absorption of higher energy electrons 
[Van AUen et al. , 1980c], or else to additional electron sources in the inner region. 

The Voyager 2 encounter produced measurements for sub-MeV electrons which app« r 
to be in conflict with the Pioneer 11 data and warrant further scrutiny. In Figure 5 the 
LECP (Low Energy Charged Particle Experiment) profiles from Krmugis and Armstrong 
[19821 are shown for 22-35 MeV electrons, >1.5-MeV electrons, and other rates for ions 
and CRAND protons in the region where the profiles in Figure 4 converge In contrast to 
the Pioneer results, the LECP data indicate an increasing low energy electron component 
in the region inwards of Enceladus. One should note, however, that the background levels 
increase proportionately to the higher energy electron rate such that electron contarmna jon 
becomes dominant inwards of Mimas. A critical question for further investigation is whether 
the relatively low background levels at 4-5 R S allow one to interpret the nominal profiles 
as being due to low energy particles. Since this question is critical to interpretation of the 
LECP data in terms of E-ring absorption effects (see accompanying section by . J, 
possibility of electron contamination is discussed further below. 

Integral Intensities and Detector Responses 

Electron intensities near a few MeV are high in the inner magnetosphere, although 
nowhere near peak levels two orders of magnitude higher at Jupiter, and there are po en i y 
valid arguments [e.g., Van AUen et al ., 1980b, c] that various different types of measurements 
at Saturn may have been adversely affected by electron contamination in the inner regmn, 
particularly inwards of Mimas. A difficulty in evaluating such arguments is that no high 
resolution measurements of the electron spectra are available in the peak flux region below 
a few MeV The best available compilation [Chenette and Stone, 1983] of electron spec 
fronT integral detector, is shown on the left side of Figure 6 for the Mimas regton, where 
contamination of the higher energy electron channels by CRAND protons was minimal, 
right side shows the measured integral intensities above 10 MeV for protons and electrons 
from Pioneer 11 and Voyager 2 near the 2.7 R S position of maximum intensity. The peak 
electron flux above 10 MeV and the peak proton flux above 600 MeV are taken from the more 
recent Pioneer 11 data analysis by Northrop and Fillius [1985], which takes into account h 
gradient-anisotropy effects relevant for the higher energy protons. 

The usefulness of data from various Pioneer and Voyager detectors is determined in par 
- by linearity of response to the intensity levels in Figure 6. Solid-state siheon detec ors op- 
erate reliably at counting rate levels up to 10 4 cps and should have had reasonably good 
efficiency for measurements of electrons above five MeV and of protons at all energies. How- 
ever any such detector directly exposed to electrons incident with lower energies near 1-3 
MeV wo^d have been heavily saturated at peak intensities, unless the sensitive detector are 
was very small (e.g., < 10 -2 cm 2 ) and optimized for magnetosphenc measurements at hig 
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Figure 4* Inbound radial profile, at Saturn for electrons above 0.040 MeV and 0 56 MeV from the Urn- 
.. Jy on Pioneer 11 [Vta Allen el .1., HWKb; Ven Allen 1984]. Note .he converge nce 
of the T two prolUe in word, from (he orbit of Tetbye, indicting that met electrons are above the higher 
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Figure 5: Inbound radial profiles from the Low Energy Charged Particle Experiment (LECP) on Voyager 
2 for low energy protons, electrons, and for MeV electrons and high energy protons. [Krimigis and Armstrong, 
1982] The solid curves indicate the sector-averaged, nominal counting rates and the small, horisontal bars 
mark the background measurements for low energy particles as the instrument moves behind a sun shield 
during its rotational movement. 
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flux levels. Typically, the detectors responding to electrons with higher energy thresholds 
were well-shielded by passive material against direct incidence by lower energy electrons. 
The unshielded detectors were either turned off, as for the Cosmic Ray Experiment (CRS) 
on Voyager 2 [Vogt et al., 1982; Schardt and McDonald, 1983], or else counted steadily at 
a constant saturation level and did not significantly affect other detector counting rates. 
Although secondary radiation processes (i.e., bremsstrahlung) may have adversely affected 
counting rates of some single detectors, particular those with minimal shielding, such pro- 
cesses probably had negligible effects on rates for multiple detector coincidences. Four-fold 
coincidences were required for ”7-17 MeV electrons” in the Pioneer-11 Chicago experiment 
[Simpson et al., 1980a, b] and for ”63-160 MeV protons” in the Voyager 2 Cosmic Ray Ex- 
periment, and electron background effects from accidental coincidences in single detectors 
were probably minimal for these rates. These two rates have interesting implications for the 
dust limit analysis as discussed further on. 

” Pathological ” Responses to Me V Electrons 

The high electron intensities at low energies make it more probable that "pathological” 
response modes may become significant in some instruments. Early concerns about electron 
background were discussed by Van Allen et al. [1980b, c] and countered to some extent 
by the responsible experimenters for the more suspect instruments [e.g., Simpson et al., 
1980b, 1981], the high-order coincidence technique being the most effective defense against 
the contamination argument with respect to electron bremssstrahlung. Another source of 
contamination is leakage of energy electrons into the region of a sensitive detector via a small 
pathway through a passive shield which should otherwise stop the highest intensity electrons. 
A recent experimental calibration at Caltech has verified leakage of electrons near and above 
three MeV through the side of the Electron Telescope (TET) in the CRS experiment [R. 
Selesnick, private communication, 1989]. The measured pulse height spectra suggest that 
internal scattering of penetrating electrons contributes a major share of the background 
leakage signal in the detector. This result will have some impact, undetermined thus far, on 
the electron microsignature analysis of Chenette and Stone [1983]. 

A analogous issue is now raised regarding potential background response in the Low 
Energy Charged Particle Experiment [LECP] on Voyager 2. In this case the magnetic de- 
flection system, used by LECP to determine allowed trajectories within the open aperture 
of the instrument for incidence of low energy electrons onto sensitive detectors (otherwise 
shielded against direct incidence), may have been rendered ineffective by internal scattering 
of the intense MeV electron flux also entering the front aperture. For both CRS and LECP, 
the usage of tungsten to stimulate secondary bremssstrahlung for high energy electron mea- 
surements (CRS) or for passive wall construction (LECP) would have greatly amplified the 
leakage effects, because of tungsten’s relatively high efficiency (near 50 percent) for backscat- 
tering of MeV electrons. Although the LECP instrument uses a sun shield in one of its eight 
rotational sectors to block the front aperture and allow measurement of penetrating back- 
ground rates, this shield may have largely attenuated the primary background component, 
electrons with peak differential fluxes near the stopping energy (1.3 MeV) of the shield, and 
produced low background readings in the data of Figure 5 for the inner magnetosphere. Due 
to its potentially adverse impact on the accompanying E-ring analysis of L. Hood, the LECP 
background response warrants further scrutiny. 

The Pioneer and Voyager measurements in the inner magnetosphere of Saturn are not 
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Ficure 7- Satellite encounter times versus energy for electrons and protons drifting longitu<Unally in the 
radial orbit IfEnceLus at equatorial pitch angles of 30 and 90 degrees [Van Allen et al., 1980c], For electrons 
peo.t*. th« ete el~ the time. fo. (e.g„ toUl) .b~P>» b„ EpeeW... 



in Figure 7. 
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sufficient to resolve the differential shape ot the Saturnian electroo spectrum m the cn . si 
enerev region near 1-2 MeV where drift resonance effects are expected. Figure 7 shows the 
cSated loss times from Van Allen et ai [1980b, c] of electrons ,n the sweeping region of 
Enceladus The large increase in loss time near 0.9-1.3 MeV corresponds to the energies at 
which electrons have minimal relative drift in magnetic longitude with 

Although the resonance would significantly distort the measured electron spectrum at 0.L3.0 
MeV, electrons at other energies and protons at all energies have finite loss times governed by 
energy-dependent velocities for longitudinal drift relative to the satellite and are unaffected. 
The flatness of the integral electron spectrum up to 1-2 MeV in Figure 6 could ^ accounted 
for by a resonance feature in the differential spectrum as predicted. Since inwardly diffusing 
electrons are accelerated, the resonance peak of the electrons getting by Enceladus moves to 
a higher energy with decreasing L towards Mimas’s orbit whereupon, peeping y 
introduces a new resonance as shown in Figure 8 from calculations by Vogt et a/. [1982J. 
Similar resonances would be produced by satellites at the orbits of Tethys an eyon . 

The pathological interpretation of the LECP data in Figure 5 is then that the differ- 
ences between the nominal counting rate levels and the background levels are due solely to 
absorption of the resonance peak electrons by the LECP sun shield during e * c 6 r ™n 
measurement. In this interpretation the nominal rates for 22-35 keV electrons and fo 
protons are dominated by internally-scattered MeV electrons inwards of E^eiadus. The 
background level also increases radially inward from Enceladus where a significan 
of the resonance electrons are above the energy required to penetrate the sun shield, and 
acceleration increases this fraction with decreasing L. Between Tethys and Ence * e 
bonafide low energy particles may be decreasing in intensity radially inward « the back- 
ground component increases, so the total nominal rates are nearly constant. The cons J 
of the background rates for 22-35 keV at 4-6 R s suggests that these rates may be dominated 

by other radiation components in this region and beyond. 

A possible ’’smoking gun” for the above pathology is found m the LECP dectron spectra 
reported by Krimigis et al. [1983]. These spectra were measured by LECP instruments fo 
b Q P th Voyager encounters with Saturn and are reproduced here as Figure 9 (their Figure ). 
In the outer magnetosphere (L > 10) both sets of spectra are either exponential or power 
law at energies up to one MeV, as one might expect for the usual magnetosphenc processes. 
In the Voyager 2 data at L = 4 - 7, however, one sees very unusual spectra with power laws 
It 30-100 keV and large humps in the spectra at 100-1000 keV. The absence of such features 
at comparable positions in the Voyager 1 data suggests a time-dependent effect, perhaps 
analogous to that seen in the Pioneer 11 electron data. 

The innermost spectrum at L = 4.2 - 5.4 shows a local maximum quite ~ 
predicted by Van Alien et oL [1980c] for drift resonant electrons diffusing past the Enceladus 
NoL that the flux falls rapidly at 1-10 MeV as also expected. In the absence of 
published spectra with better spatial resolution we can only speculate that the maximum 
rows larger inwards of Tethys’s orbit (4.9 R s ). The other Voyager 2 spectrum for L - 6.9 - 
5.9 Rs does not show the same peak structure but may be partially affected by a Dione 

16 3 RqI drift resonance. , . . , . 

' Although no similar LECP data are given inwards of Enceladus, probably due to un- 
certainties in background corrections, we further speculate that the resonance P eaks “ 
electron spectra become totally dominant there, and that the low energy fluxes cannot be 
measured without quantitative corrections for background from electrons m these peaks. In 
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the intermediate region at 4 - 5 R S) where local minima appear in the radial profiles for 
low energy electrons and ions, the electron spectra are probably an even mixture of the 
soft, low energy component and the high energy, resonance component. The apparent in- 
ward increase from the transition region of the low energy counting rates may simply reflect 
increasing background levels from the now dominant resonance component. 

Differential Spectra for Protons 

Differential flux measurements of proton energy spectra could provide a sensitive measure 
of absorption in diffuse dust layers near the ring plane, since the proton loss rate in dust 
increases with decreasing energy below 300 MeV due to ionization energy loss. In Figure 1 
we compare the calculated CRAND source spectrum for protons from Cooper [1983] with the 
very limited data on spectral measurements from Pioneer 11 and Voyager 2. Relative to the 
nearly flat source spectrum below 100 MeV, the measured spectra from the Voyager 2 Low 
Energy Charged Particle (LECP) experiment [Krimigis and Armstrong, 1982] show strong 
depletion of low energy protons at L = 2.73, possibly indicating dust-dominated ionization 
losses in that region. A comparable differential flux at 63-160 MeV is also available from the 
Voyager 2 Cosmic Ray (CRS) experiment of Vogt et al. [1982]. CRS proton measurements at 
lower energies (48-63 MeV) show apparent increases of differential fluxes from higher energies, 
but spectrum-dependent corrections for electron background make this determination highly 
uncertain [Schardt and McDonald, 1983]. We assume that the LECP data better represent 
the actual shape of the spectra below 100 MeV. 

Overall, the differential fluxes in Figure 10 are consistent with a CRAND source spectrum 
and energy-dependent losses below 100 MeV from interactions with dust. The LECP spectra 
at 10-100 MeV are a reasonable match to a theoretically derived spectrum [Schardt and 
McDonald, 1983] of the form J(E) « E 1+7 , where 7 « 0 approximates the source spectrum 
below 100 MeV. At higher energies, from 80 MeV to above 600 MeV, the integral flux data 
from Pioneer 11 [Fillius et al., 1980; Northrop and Fillius, 1985; Van Allen et al., 1980b] 
are consistent with the indicated differential spectrum (7 ~ -3) in Figure 5 and with the 
CRAND source spectrum. 

Since radial diffusion broadens the effective region for observable effects on the proton 
spectra, it is plausible that the Voyager 2 spectra at 2.73 Rs were affected by absorption 
occurring further outwards in the G ring. However, the Voyager 2 LECP found similar 
spectra at 3.4 R s , although with less depletion near 20 MeV, so one cannot rule out local 
absorption effects from the measured spectra which actually appear quite consistent with 
expected losses in dust. One might speculate that the dust densities are just high enough 
to affect the lower energy protons, but that the integral intensities of higher energy protons 
are minimally affected. Since our CRAND modeling of radial profiles for the latter shows no 
significant dust effect, the marginal column densities for lower energy losses would correspond 
to the ” worst case” limits for 100-MeV protons. These limits are relative high (10 gm/cm ) 
in the E ring, and the same limits would give observable effects on the radial profiles near 2.7 
R s in contrast to the modeling results. Alternatively, some other explanation may be found 
for the spectra which is unrelated to dust effects (e.g., low energy anomaly in the CRAND 
source spectrum, unaccounted-for effects of background contamination?). 

Pitch Angle Anisotropies 

Anisotropy measurements for the pitch angle distributions of trapped protons are also 
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ENERGY (MeV) 


Fi eure 10: Calculated [Cooper, 1983] and measured (dark circles from Krimigis and Armstrong [1982]; 
open circle from Vogt et «/., 1982) differential energy spectra for CRAND protons at L « 2.7. The higher 
energy dthed spectrum is ^rapolated from the integral fluxes above 80 MeV and 600 MeV » Figure 6. 
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sensitive to the presence of diffuse dust or larger aggregates of co-orbiting material near the 
ring plane. Voyager dust impact data [Scarf et al., 1982, 1983] indicate that detectable dust 
layers at Saturn are confined within a few thousand kilometers of the ring plane. Protons 
with large equatorial pitch angles ao at the ring plane have larger pathlengths through a 
layer of uniform thickness and therefore have greater losses. The absorption of the equatori- 
ally trapped proton population becomes especially prominant when observed in the CRAND 
proton population. Such protons are nominally expected [Blake et al., 1983] to have prefer- 
entially higher injection rates for large values of ao. 

In principle, the correlation of anisotropies and differential energy spectra should allow 
differentiation between absorption effects due to dust and to other absorbers, including major 
satellites. Declining LECP spectra below 100 MeV and reduced CRS anisotropies are jointly 
indicative of absorption in a layer of small dust particles, where the spectral effect arises 
from slow energy losses as compared to catastrophic single impacts with large particles (i.e., 
radii much larger than ten centimeters). The random distribution of larger particles (e.g., 
moonlets up to kilometers in size) over a wide radial range in a vertically thin layer can also 
produce anisotropies like those observed, but only the smaller particles produce the spectral 
effect. 

Figure 11 shows a Voyager 2 measurement of proton anisotropies for for 63-160 MeV 
protons as reported by Schardt and McDonald [1983]. Ratios of inbound and outbound 
fluxes at different magnetic latitudes were used to determine the ’’pancake” distribution 
parameter n for the differential distribution j(ao) oc sin n ao of directional flux near the ring 
plane. The largest values of n (i.e., more anisotropic fluxes) were found in the highest proton 
flux regions, as expected for nominal geometry of the CRAND source [Blake et al., 1983], 
while low values (i.e., more isotropic fluxes) were found in the vicinity of the G-ring, Mimas’s 
orbit, and in the E-ring region beyond 3.3 Rs- In the latter region, the radially outward 
decrease in n is suggestive of a strong E-ring effect but may have other explanations. 

One such explanation was first advanced by Cooper [1983], who noted that re-absorption 
of CRAND-parent neutrons would occur preferentially for emission from ring source regions 
at more tangential angles with respect to the ring plane. Figure 12 shows Cooper’s calcu- 
lations for the dependence of CRAND injection efficiency (i.e., the injection coefficient) on 
equatorial pitch angle with a simple, thick slab model the parent neutron source in the main 
rings. Although this calculation underpredicts the anisotropy of the injected protons inwards 
from 3.3 Rs, where the discrete nature of the neutron-emitting bodies in the rings [Blake et 
al., 1983] accounts for the largest anisotropies observed, it does predict decreasing injection 
efficiency at large distances for protons with large equatorial pitch angles. This effect arises 
because (1) the parent neutrons have increasing slant pathlengths for escape from the source 
slab, and (2) the radial distance from the source to the injection point increases faster for 
injection near the equatorial plane at large L than at high latitudes, where injection sites 
remain directly over the rings for the smallest equatorial pitch angles. Whether this to- 
tally explains the isotropization of protons in the E-ring region is unclear, but more realistic 
modeling of the ring source would be helpful to resolve this question. What is surprising is 
that no one, myself included, has addressed this rather obvious problem since Schardt and 
McDonald reported the anisotropy results in 1983. 

The pitch angle distributions of injected protons should be sensitive to size distributions 
and vertical positions of neutron-emitting bodies within the main rings, where the typical 
sizes of such bodies are of order one meter in radius as compared to 10-100 meters for the ver- 
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Figure 11: Anisotropy index n for 63-160 MeV protons from the Cosmic Ray Experiment (CRS) on 
Voyager 2 [Schardt and McDonald, 1983]. Large values indicate relatively more eqnatorially trapped protons, 
while smaller values arise for more isotropic distributions. 
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tical ring thicknesses [Tyler et a/., 1982]. At some radial distance, dependent on the ratio of 
interparticle separation distances and the vertical dispersion scale, one should «P«tstrong 
reabsorption of tangentially scattered neutrons, perhaps leading to the observed effec 
yond 3.3 R s . At higher proton energies, above a few hundred MeV, this effect should become 
more pronounced, since scattering angles decrease with increasing secondary neutron energy. 

In this regard, improved measurements of CRAND proton energy spectra and anisotropies 
in this rein by Cassini might therefore used to probe the internal structure of ring particle 
distributions in rough analogy to an nuclear scattering experiment m high energy physics. 
Sophisticated monte carlo codes, originally developed for work in high energy p ysics, can 
be used to explore these possibilities. 

High Energy Electrons ,, , , 

A unique signature of diffuse dust in regions of peak proton intensity would be the 
presence of highly relativistic electrons with radial intensity profiles roughly approximating 
those of the parent CRAND protons. Since longitudinal drift rates of both protons and 
electrons scale roughly with kinetic energy, one expects very low flux levels for electrons 
above 10 MeV as observed [Chenette and Stone, 1983] in satellite regions where e pro on 
fluxes are low. The electron and proton fluxes may be much higher ebewhere m the absence o 
strong absorption, but only if local in-situ sources are active. Since CRAND already provi es 
the required source for the protons, nuclear interactions of the protons above 300 MeV would 
provide a source of high energy electrons via pion production at a level proportional to the 
proton fluxes and to the column density of diffuse dust. Upon injection, the pionic electrons 
would have energy spectra extending into the 10-100 MeV range, the mean pion-muon ecay 
energy for electrons being near 50 MeV. If the trapping lifetimes of these electrons were 
limited largely by radial diffusion, as for the protons, the observable spectra of the time- 
averaged electron population could be comparable in shape to the injection spectra, while 
a more significant loss rate from dust interactions would tend to soften the time-averaged 
spectra due to the increasing losses from bremsstrahlung at higher energies. In either case 
one would expect significantly flatter spectra above 10 MeV from the dust source than are 
observed in satellite regions, where inwardly diffusing electrons would dominate the integral 
fluxes below 5-10 MeV, and the local source component would be rapidly reabsorbed. 

Experiments which do not discriminate between penetrating P^tons and the higher 
energy electrons can produce highly ambiguous results as seen m the 7-17 MeV Elec * ron 
profiles in Figure 13 from the University of Chicago’s . ringed parhcle eqemwfi .onboard 
Pioneer 11 [Simpson et al, 1980a, b, 1981; R. B. McKibben, unpublished data, 1988]. These 
profiles may be considered to be a mixture of contributions from electrons at energies above 
7 MeV and from protons above 30 MeV. The electrons are dominant in radial zones of proton 
absorption by Janus/Epimetheus near 2.5 R S and in the vicinity of the F-nng region, while 
the protons account at least partly for the peak counting rates near 2.4 R s and 2.7 R s . 
The large change in counting rates in the vicinity of the absorption regions is probably due 
to rapid changes in the electron intensities over the intervening four-hour interval between 
inbound and outbound observations [Fillius et al ., 1988]. 

The outbound drop in electron intensity over a wide radial range in Figure ^reveals 
the rising shoulders of profiles for more stably trapped particles winch are strongly absorbed 

within the satellite and ring regions. An outstanding question .. the “ 

stably trapped particles in the peak intensity regions. The possibility remains that electrons 
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Figure 12: Calculated injection coefficients versus equatorial pitch angle for CRAND protons at selected 
radial positions [Cooper, 1983]. These coefficients were calculated for parent neutrons emitted from sources 
ii a .r,A tt rinn modeled as a uniformly thick annular disk at 1.5-2.3 Rs. Note the reduced injection 


in the A and B rings, modeled as a uniformly thi 
at large pitch angles near Enceladus and beyond. 
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RADIAL DISTANCE IN SATURN RADII (Rs - 60,000 km) 


Figure 14: Inbound (top) and outbound (bottom) radial profiles from Pmneer 1 for petrous above 
10 MeV and protons above 600 MeV, as determined by the gradient-anuotropyanalysisof North* p 
FUlins [1985] Note the significant inbound-outbound asymmetry in the electron flux profiles inwards of 2. 
R S which is similar to that observed for the "7-17 MeV electrons" in Figure 13. 




at 10-20 MeV dominate the counting rates in these regions. In this case the peak intensity 
regions observed outbound might be partislly, if not totally, dominated by higher energy 
secondary electrons from the proton-dust interactions. Although the pionic electrons would 
be produced with average energies near 50 MeV from pion-muon decays, long term trapping 
and radiation losses in the dust might produce an equilibrium population of electrons with 
lower average energies in the 10-30 MeV range. Thus this secondary population might 
contribute substantially to the 7-17 MeV electron counting rates. On the other hand, we 
cannot yet rule out the possibility that this rate is entirely due to protons without new 
measurements by Cassini. 

Although the high energies of the pionic electrons and their proton parents make dis- 
crimination between the electrons and protons difficult, the gradient anisotropy effect has 
already been used in the Pioneer 11 data analysis [Fillius and McDwain, 1980; Northrop 
and Fillius, 1985] to estimate this separation and to establish the presence of electrons at 
energies above 10 MeV with radial intensity structure similar to that of protons above 600 
MeV. The electron fluxes are separated from those of protons by the negligible east-west 
anisotropies of the electrons in comparison to the very significi it anisotropies observed for 
protons above 600 MeV. The latter have large gyroradii which offset guiding centers by 10 3 
km for protons measured in regions of large radial intensity gradients. Figure 14 shows the 
deconvolved profiles from Northrop and Fillius [1935] for the electron and proton components 
as measured (a) inbound and (b) outbound by Pioneer 11. Although the electron intensities 
in the 2.5 Rs region of satellite absorption show a large outbound drop in intensity similar to 
that of the ”7-17 MeV Electrons” in Figure 13, the intensity peaks of both the electrons and 
the protons change little in magnitude or position between inbound and outbound observa- 
tions and are probably associated with stably trapped populations. Although proton-dust 
interactions were not considered by Northrop and Fillius in their discussions of potential 
electron sources, it entirely possible that the electrons near 2.7 Rs are pionic in origin and 
therefore are also strongly indicative of a radially extended diffuse dust layer not unlike the 
more easily observable E-ring material beyond the orbit of Mimas. 

Microsignatures 

Among the most controversial aspects of the Pioneer and Voyager encounters were un- 
doubtably the simultaneous detections of microsignatures by different experiments on the 
same spacecraft in the orbital region (also called the macrosignature region) of Mimas. Van 
Allen et al. [1980c] argued that the single Pioneer microsignature for Mimas (observed in- 
bound) could be explained by the effect of MeV electron response in all affected data channels 
to a direct drift shadow from Mimas, located 56 degrees eastward from Pioneer 11 in mag- 
netic longitude (MeV electrons drift westward). In rebuttal, Simpson et al. [1980b] argued 
strongly that the observed microsignature, evident in Figure 15 from counting rates nomi- 
nally for 0.5- 1.8 MeV protons and 7-17 MeV electrons, could not be explained by electron 
contamination. 

The so-called ”L1NL2” counting rate for the MeV protons shows highly aliased response 
in some locales within the Saturnian magnetosphere. The high intensity peaks, evident in 
Figure 3 for the nominal MeV (L1NL2) proton profile, have been attributed [Simpson et al., 
1981; Cooper, 1983] to side-penetration of the associated detector by CRAND protons at 
much higher energies, while the low but finite counting rates in the Mimas region are still 
attributed to MeV protons, based solely on the proton-like nature of the measured energy 
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SPACECRAFT TIME (UT). DAY 244. 1978 S481 

Figure 15: The "Mimas" microsignature as observed inbound in the "7-17 MeV electron” and 0.5- 
1.8 MeV proton” counting rates from the University of Chicago experiments on Pioneer 11 [Simpson et 
al., 1980a, b). Increases in both rates on each side of the Mim as sweeping region are most likely due to 
background from high energy protons. This feature was also detected in the Pioneer 11 data of Van Allen 
et al. [1980a, b,c]. No corresponding feature was detected outbound. 


Time of day 238,1982 



Figure 16: Inbound (left) and outbound (right) counting rates in the Mimas region from the Voyager 2 
Cosmic Ray Experiment (CRS) [Vogt et al., 1982; Chenette and Stone, 1983) The "BS4” rate is a double 
coincidence for detectors sensitive to 2.5 MeV and 6.5 MeV electrons. The ”5 MeV electron rate is also 
sensitive to electrons near three MeV due to side penetration response in the CRS Electron Telescope [R. 
Seles nick, private communication, 1989). Spacecraft maneuvers may have affected the inbound profiles, while 
the "Mimas ghost” microsignature is observed only outbound near the outer radial limit of Mimas s eccentric 
orbit, which appears to sharply define an absorption boundary for the high energy protons. 
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deposits in the detectors. On the other hand, it is clear that the L1NL2 rate will count 
electrons at some level, since this rate is correlated to the inbound-outbound flux decrease 
for electrons in the 2.5 Rs region, although no pulse height data are available from there. 
There was no appreciable decrea^t m the Mimas region. Ongoing analysis of calibration 
data from laboratory tests of a similar detector with a high flux, beta electron source [J. F. 
Cooper and J. P. Wefel, unpublished data, 1989] may produce a more definitive evaluation 
of the electron background problem for this detector. 

Voyager 2 detected a single, outbound microsignature, shown in Figure 16, m the Mi- 
mas region in both electron [Vogt et a/., 1982] and low-energy ion [Carbary et ai, 1983] 
channels, this time when Mimas was 147 degrees west of the observing spacecraft. This was 
interpreted by Chenette and Stone [1983] as supporting the original view of Simpson et ai 
[1980b]. Chenette and Stone further proposed that the observed microsignatures probably 
arose from near encounters with 100-km clumps of co-orbiting material which may have been 
widely distributed throughout the macrosignature region to give the observed, fifty- percent 
probability of detection as microsignatures. Carbary et al. also found a single inbound mi- 
crosignature at the orbit of Enceladus in low energy ion and MeV electron data. They also 
attribute the Enceladus feature to clumping of co-orbiting material, perhaps associated with 
the E-ring. The evidence for clumping is even more dramatic in the F-ring, where Cuzzi and 
Burns [1988] have argued that three of the five F-ring microsignatures, observed by Pioneer 
11 experiments [Simpson et al . , 1980a, b; Van Allen et al L, 1980a, b] and shown in Figure 
17, could not have been produced by the F-ring itself and must have arisen from random 
encounters with clouds of debris in the general vicinity. 

Interpretation of the microsignature observations may lead to new planetary scierce 
but also indicates potentially grave hazards for the Cassini orbiter in the satellite regions. 
The common "horseshoe” orbits of the fragmented satellites Janus and Epimetheus at 2.52 
R s are the most easily observed example of stable orbital configurations in which many 
bodies can occupy the same radial zones by means of N-body gravitational interactions [e.g., 
Cuzzi and Burns, 1988; J. Burns, private communication, 1989]. Assuming that no residual 
problems remain with the microsignature observations, one is forced to conclude that large 
numbers (e.g., thousands) of localized debris clumps may co-orbit with the F-ring and the 
larger satellites. On the other hand, some of the microsignature data may be subject to 
"pathological” responses to high intensity background from MeV electrons, viz-a-viz the 
concerns of Van Allen et al. [1980b, c], so further experimental and analytical effort would 
be appropriate to evaluate the microsignature data and the likelihood of such responses. 

3. CRAND Models for the Proton Observations 

The thesis work of Cooper [1983] provides our starting point for modeling the CRAND 
proton observations in the inner magnetosphere of Saturn. This work modeled ^the main 
rings as a solid slab of material with an average column density of 100 gm/cm , a value 
in reasonable agreement with Voyager measurements from microwave-radio [Tyler et ai, 
1982; Marouf et ai, 1983] and starlight [e.g., Esposito et ai, 1983] occultation data. This 
slab was modeled with monte carlo calculations as a target layer for interactions by high 
energy cosmic rays which can penetrate Saturn’s magnetosphere and strike the rings at 
energies above 10-20 GeV/nucleon. The resultant interactions produce secondary radiation 
including neutrons, which then become parents via beta decay for the trapped high energy 
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Figure 17: Inbound (A) and outbound (B) counting rates from Pioneer 11 [Simpson et al . , 1980a, b] 
for MeV electrons in the F-ring region where five microsignatures are evident. The flatter shapes of the 
outbound profiles are associated with the general drop in electron intensity inwards of 2.6 Rs- The "D5N6 
profiles correspond nominally to "7-17 MeV electrons” from the University of Chicago’s Main Telescope, 
while the other two profiles are accumulated tom low-order coincidences of lower energy electrons in the 
front detectors of the same instrument. 
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protons. An important constraint on Cooper’s original model was supplied by . irect j? 

H observations of secondary charged particles and gamma rays from these mteractions m 
the rings [Chenette et ai, 1980; Cooper et ai, 1985; Van Allen and Randall, 1985]. These 
observations estabhshed that the observed secondary radiation intensities were consistent 
with results from the slab model and that the neutron enussion could be accurately modeled. 

Since free neutrons have average decay lifetimes of 10 seconds, only a very sm 
decay within the magnetosphere, so the injection rate of trapped, decay protons into the 
radiation belt regions is very small. We will use calculated profiles in Figure 18 of neutron- 
d«.y proton injection from Cooper-, work to deSne the rndinl .ource profiles shown o 
different proton energies. The nearly exponential shapes of these profiles up to several 
hundred MeV allow simple parametric forms to be used for radial dependence of the sou 
rates, which vary by factors of two or less between adjacent trapping regions. For narrow 
radial zones a constant source rate can be assumed. 

L03S ()nce injected, the trapped protons undergo radial diffusion at a rate determined by 
model fits to observations. Large trapping times (* 10 - 10 years) are required to bull up 
significant fluxes at the levels observed. In the model the sweeping by satellites is treated by 
the definition of totally absorbing boundaries at the approximate positions of the satellites 
where allowances are made in the fitting procedures for the finite gyroradu of the proton 
which increase the effective cross-sectional areas of the satellites ^r sweeping. Simple analytic 
forms for the solutions of radial transport equations with the CRAND source term, ra i 
diffusion, and the satellite boundaries for each diffusion region are obtained in the hmit tha 
the proton profiles are time-independent. A principal limitation of the original modelis that 
no attempt is made to include the obvious dust-related losses in the region of the G ring, 
the only modeled losses being those due to radial transport into regions of strong absorption 
from satellites and or rings, where the proton flux can be assumed to approach zero. 

Figure 19 shows the results from Cooper [1983] for fits to University of Chicago profiles for 
high energy protons in the inner magnetosphere of Saturn. The energy-dependent response 
functions [Simpson et ai 1981; R. B. McKibben, unpublished data, 1982] of de ) cctors 
the observed profiles have been folded into the model parameter fits. In general the diffusive 
gradients of the proton profiles are well-matched by the model curves in the vicinity of he 

satellite boundaries for a diffusion coefficient of the form Dll - A)L > wbere n r ° m 
fit and n = 10 for diffusion driven by magnetic impulses from magnetosphenc disturbances. 
In the region around the proton peak at 2.4 R S the diffusive profile matches the observed 
reasonably well, while the fit is clearly not satisfactory m the G-nng region where dust b ses 
have not been included in the model. Beyond 3.3 R S the excess in the calculated rate relative 
to that measured indicates either that the source function requires correction or that E-n g 
losses should be included in the model. In principle, the profiles and relative “tensities 
the three major protons peaks can be totally explained m terms of the radial width of the 
lossless regions around the peaks and by the radial dependence of the source > n J ect, °" 
and of the diffusion coefficient. On the other hand, the observational and model-rela 
uncertainties dictate that we estimate upper limits for potential losses having no o vious 
effect on the observed radial profiles. 
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Figure 18: Radial profiles from Cooper [1983] for integral, omnidirectional injection rates of CRAND 
protons above selected energy thresholds. 



Figure 19: Best fits of the Cooper [1983] CRAND model to Pioneer 11 radial profiles from the University 
r rh * f or protons above 30 MeV. The fitted curves correspond to models including truncation of the 

nnd d.tWO, functions Wo. « M«V (1, n.d »0 M.V (2), - - 

latter was found to be most consistent with observed ratios of the two detector responses. The mode c 
* fitted to thc observed absorption boundaries and the maximum proton intensities. Model fits to the 

* peaks are conSZt with a diffusion coefficient D, L m 1Q-«X« Rs’/‘« for P-tons > 10* MeVm 
three pe n( . tos _ here The Cooper model did not include effects of proton absorption in the E nng. The 

» a., *w <« * ■ *, d-*- . m .» 

changes in the source and injection geometry. 


28 





Radial Transport with CRAND and Losses 

Using the CRAND model to set limits on dust, abundances requires modification of the 
Cooper model to allow for distributed losses througout the diffusion region. The effects 
of satellite sweeping can still be retained in the boundary conditions as before. The most 
tractable solutions for the modified model arise when the losses are approximated by an 
average time scale for total loss which may be an analytic function of radial position and 
proton energy. Such solutions have previously been used [e.g., Hood, 1983, 1985] for lossy 
diffusion problems not including distributed sources. Work done in collaboration with Lon 
Hood has now generated the required solutions for distributed sources and losses. Develop- 
ment of these new solutions allows estimates to be made on upper limits for dust abundances, 
based on model fits to the observed radial profiles for high energy protons. Special cases 
of these solutions for the problems at hand are described in Hood’s accompanying report. 
More detailed versions will be published later by Hood and myself. 

The radial transport of charged particles in any stable, dipolar magnetic field is described 
by the following differential equation: 



This equation [e.g., Schulz and Lanzerotti, 1974] governs the time and spatial evolution of 
the phase space density function f(M,K,L) with respect to the three adiabatic invariants of 
charged particle motion: the magnetic moment M for gyromotion, the latitudinal bounce 
invariant K, and the radial parameter L defining the local field line. For simple calculations 
the usual procedure is to set K to zero in the limit of equatorially trapped particles, while 
allowing M and L to vary freely. Since the trapping times of CRAND protons are very long, 
the time dependent term df/dt is assumed to be negligible. The phase space space density 
is related to the differential particle flux j by f = j/p 2 , where p is the particle momentum on 
the local field line. 

In the case of distributed injection by CRAND, there is little net acceleration, since 
protons diffuse both inwards and outwards from the intensity peaks, so we define an approx- 
imately constant momentum. In this case the phase space density is roughly proportional to 
the local integral flux at the level of the present analysis. It may also be more appropriate 
in such a case to use a one-dimensional transport equation of the form 



which can be easily solved in the simplest cases for constant values of Dll, S, and r - These 
solutions are discussed in Hood’s report. 

The last three terms in the transport equation pertain to radial transport with diffusion 
coefficient Dll, the CRAND source injection rate 5, and the loss rate L due to dust interac- 
tions and/or other processes. As for the phase space density, the source and loss terms are 
proportional (when multiplied by p 2 ) to the source and loss rates in flux units. We make 
further simplications by defining constant values for the source and loss terms, reasonably 
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The actual dust abundances could be quite close to the upper limit values and still 
have only marginally observable effects on the radial profiles for integral counting rates of 
CRAND protons, due to the competing losses from radial diffusion which steepen radial 
gradients in shoulder regions adjacent to intensity peaks. Hood’s quantitative fits sdiow 
that the ’’nominal” dust limits are a factor of three (i.e., an e-folding factor m loss time 
units) below the worst case values when the radial profiles are fitted with a one-dimensional, 
lossy diffusion model. On the other hand, these first order fits do not include effects of 
radial dependence in source and diffusion parameters, radial variations m proton energy and 
pitch angle distributions, or the effect of finite gyroradii (e.g., the gradient-anisotropy effect 
analysed by Northrop and Fillius [1985]). Thus, the worst cases values in Table 1 should be 
regarded as the more conservative estimates for impact risk assessments. 

In comparison, the very substantial proton modulation by the G-ring (c.f., Figure 3) is 
associated with an estimated G-ring column density of 2x10 7 gm/cm , based on Hood s 
update of earlier modeling work done for the G ring by Van Allen [1983, 1987]. If the dus 
column density near the proton peak at 2.7 Rs is near the worst case value m Table 1, the 
higher proton intensity near the peak may actually lead to rates for proton interactions in 
dust which are within a factor of three of those in the G-ring region. Radial diffusion may 
also broaden the radial extent of regions effected by absorption in adjacent regions of stronger 
absorption. In this case, observable effects of absorption would be evident in more sensitive 
measurements over a broad radial range, even outside the principle absorption regions. 

The Satellite Regions 

Evaluated quantitatively, the microsignature observations yield very high estimates for 
column densities of co-orbiting material in the regions of observation. For the Mimas ghost 
microsignature we get a typical column density of order lxlO" 2 gm/cm for each clump of 
diameter 10 2 km by using parameters calculated by Chenette and Stone [1983]. Average 
over the width of the Mimas region (0.1 R S ) and the number of clumps (« 3x10 ) consistent 
with the observations, the ”macrosignature” column density becomes 4x10 gm/cm . Al- 
though microsignatures in the Janus/Epimetheus and Tethys regions were clearly related to 
near encounters with known satellites (e.g., Van Allen, 1984), the observational statistics are 
such that we cannot rule out comparable amounts of co-orbiting matenal m re S lon *' 

Another quantitative approach was first discussed by Schardt and McDonald [1983], who 
found that the measured fluxes and anisotropies for 63-160 MeV protons in the Mimas re- 
gion could be consistent with absorption by Mimas alone. Their emulations * he 

reasonably comparable CRAND sources from Cooper [1983] and Blake et al [1983], the- 
oretical derivations for proton sweeping by the satellite, and further assumed that proton 
lifetimes in this region were totally dominated by satellite losses and were unaffected by 
radial diffusion. Although they reported a four-fold coincidence measurement for the mini- 
mum proton intensity in the Mimas region, it may be assumed that there is some electron 
contamination in this measurement, the electron intensities being many orders of magnitude 
higher, so there remains a substantial margin for effects by other co-orbitmg matenal. Our 
own analysis gives ” worst case” lower limits of 10" 6 to 1(T 5 gm/cm 2 for co-orbiting ma enal 
absorbing CRAND protons at 2.35, 2.5, and 3.1 R S and are consistent with the dust clump 

interpretation. 
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justifiable, since the source rate varies slowly with L, and we can only determine limits on 
the column density of a uniformly-thick dust layer in each diffusion region. The inner and 
outer edges of the strong absorption regions (i.e., orbits of principal satellites and the F ring) 
define zero-flux boundaries for the proton trapping regions (i.e,. proton gaps in our cuFrent 
parlance viz-a-viz Cassini aimpoints). 

If one defines a characteristic time scale r for dust losses, the loss term becomes L={/t. 
In the limit of lossless diffusion, r is set to infinity, while it assumes finite values for lossy 
diffusion. For proton losses in a vertically thin dust layer near the ring plane, the loss time 
scale is approximately 


R^Tgcosa o 


■where a is the dust column density in units of gm/ cm 2 , R« is the proton stopping range in the 
same units, Tg is the latitudinal bounce period in Saturn’s dipole field [e.g., Thomsen and 
Van Allen, 1980], and the cosine factor ss 0.5 accounts for the inclined pathlength through 
the dust layer. We use stopping ranges in water ice of 7.7 gm/cm 2 for 100 MeV protons and 
50 gm/cm 2 for 300 MeV protons. 

4. Worst Case Limits for Dust Column Densities 
High Intensity Regions 

Our "worst case” limits for column densities of dust are determined simply by setting the 
integral source rate 5J p /5t equal to the probabilistic loss rate J p /t, where J p is the locally 
measured integral flux of protons at any given position in L, and then solving for <t in terms 
of the other, known parameters. This limit corresponds to the rather extreme assumption 
that the proton lifetimes are totally governed by losses in dust. Since it also corresponds 
to the critical values of <r at which the integral intensity extremes are either maximal (5 
> L) or minimal (5 < L), the computed values of a are upper limits near intensity peaks 
and lower limits near intensity minima. The computed limits, upper or lower, are shown 
in Table 1 for selected L positions. The proton fluxes (J p ) were computed from Pioneer 11 
measurements of Fillius ct al. [1980] for protons sampled near the magnetic equator (i.e., 
K ss 0) at energies above 80 MeV [c.f., Figure 3]. Since the form of the differential spectrum 
is not well-determined at 100-300 MeV, we show calculated limits for 100 MeV and 300 
MeV protons, assuming that one or the other is characteristic of protons above the integral 
threshold energy. 

The inverse proportionality between the calculated <r limits and the measured flux pro- 
duces the lowest upper limit of 2-3xl0 -8 gm/cm 2 at the position of highest measured proton 
flux, L « 2.65. The other two intensity peaks, located at 2.42 Rs and 3.3 Rs, give compa- 
rable upper limits of l-2xl0 -7 gm/cm 2 , an order of magnitude less than those at the most 
intense peak. If one assumes a uniform dust layer throughout the entire inner magneto- 
sphere, except perhaps in the satellite regions where the calculated lower limits are 10 1 to 
10 2 times higher, one might take the 2.65-R s value of cr as the best-determined upper limit, 
simply because J p is highest at that position. By this reasoning, the probable effects of dust 
in the two lower intensity peaks should therefore greatly reduced relative to the potentially 
marginal effect at the highest peak. 
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Table 1 

CRAND Proton Limits for Dust Column Densities 


L 

(t (E=100) 
[gm/cm 2 ] 

<r (E=300) 
[gm/cm 2 ] 


2.35 

> 3.5x10"® 

> 5.3xl0“ 5 

(F ring) 

2.42 

< 9.6x10"® 

< 1.3xl0" 7 

(Proton gap) 

2.5 

> 2.8x10"® 

> 3.7xl0" 5 

( J anus /Epimetheus) 

2.65 

< 1.8xl0~ 8 

< 2.7x10"® 

(Proton gap) 

2.80 

> 5.7x10"® 

> 7.5x10"® 

(G ring) 

3.1 

> 2.8xl0 -6 

> 3.4x10"® 

(Mimas) 

3.3 

< 1.6xl0" 7 

< 1.8xl0" 7 

(Proton gap) 

3.5 

< 3J5xl0" 7 

< 3.9xl0" 7 

(E ring) 

3.7 

< 1.3x10"® 

< 1.5x10"® 

(E ring) 


5. Justification for Relative Merits of Cassini Aimpoints 
Proton Intensities 

Our worst case limits for column density of dust are strongly correlated (inversely) to the 
local proton intensity. The only other well known parameter is the CRAND proton source 
which varies only by a factor of two between the two aimpomts. Since the measured proton 
intensity is a factor of ten higher at 2.7 R s than at 3.2 R s , I would feel safer going in a 
the inner point, just because the higher proton intensity there is a better probe of the dust 
density. The same argument applies to a comparison between 2.7 R s and the innermost 
"proton gap” centered on 2.42 Rs, the first being preferable to the second. 

In general, the relative intensity levels in all three gaps are consistent with our expecta- 
tions for effects of CRAND and lossless diffusion as demonstrated in Figure 19. The increase 
in the radial diffusion rate by a factor of five (for D LL oc L 9 ) at 3.2 Rs gives a relative^ 
stronger diffusive loss effect which mostly accounts for the lower intensity An additional 
factor of two from radial dependence of the source rate further accounts for the order of 
magnitude difference in intensity between the peaks at 2.7 and 3.2 R s . Similarly, the lower 
diffusion rate, and the slightly higher source rate, in the 2 42-R S gap compensate for the 
relatively smaller radial width of this gap and account for the observed intensity. 

About all the integral intensity profiles can tell us is the upper limit on a uniformly thick 
dust layer, except perhaps in the satellite regions, throughout the entiremner magnetosphere. 
The best constraint on this is determined at 2.7 R s where a <2-3x10 gm/cm for th 
300 MeV protons dominating the measured integral intensities. If the dust effect there is 
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marginal (i.e., within a factor of three of the worst case limit) , there should be less observable 
effects in integral intensity profiles in the other proton gaps. For these the lower intensities 
give correspondingly lower proton impact rates on dust relative to loss rates from radial 
transport into the strong absorption regions of the satellite orbits. 


Proton Energy Spectra 

Although the CRAND model fits to integral intensity profiles are consistent with lossless 
(i.e., no dust) diffusion at 2.7 R s and elsewhere, higher sensitivity might be achieved by 
looking at protons below 100 MeV where ionization energy losses are more signify nt. In- 
deed, some data from Voyager [Armstrong and Krimigis, 1982] indicate dropoff towards Me 
energies in spectra measured near 2.73 R S and near 3.4 R s . Other spectra from Schardt and 
McDonald [1983] do not show this dropoff, but probably have significant electron contamina- 
tion at the lower energies (27-63 MeV as compared to 63-160 MeV). From these data I infer 
a significant dust absorption effect for the lower energy protons. Since the proton energy 
loss range increases by an order of magnitude from the lower energies to 300 MeV, the dust 
effect might be significant for the lower energies while still being very marginal at the higher 
energies which dominate the integral counting rates used for the radial profile fits. Thus 
I doubt that the real dust column density is much less than an order of magnitude below 
the worst case limit for 2.7 R s , when the spectral information is considered. Furthermore, 
the potential evidence for secondary electron production in the proton gaps argues for some 
finite, but yet to be determined amount of diffuse material. On the other hand, there may 
other explanations for the proton spectra and high energy electron results, so these cannot 
be considered as very definitive at this point in terms of providing information on dust bmits. 


Proton Anisotropies 

The anisotropy peak at 3.2-3.4 Rs in the Schardt and McDonald [1983] data certainly 
suggests that this region is relatively clear of absorbers as compared to the G-nng and the 
Mimas region. Since one is still apparently within the influence of the G-nng via diffusion 
at Voyager 2’s innermost L-shell of 2.73 R s , the Voyager 2 measurement there may not tell 
us anything about the actual anisotropy near the peak intensity level at 2.65 R s . However, 
Van Allen ei al. [1980b] did note that roll modulation in their Pioneer 11 counting rates 
was consistent with a high degree of anisotropy (n w 4), but they did not specify the radial 
dependence. Close examination of Figure 4 in their Science paper [Van Allen et al , 1980a] 
shows no significant variation in the roll modulation of the rates for >80-MeV protons in 
the region at 2.6-2.8 Rs, so it is unclear whether the anisotropy increases from a low level 
at 2.8 R S in the Voyager measurement (c.f., Figure 11) to the high level estimated from the 
Pioneer data further inwards. One might also note that the roll modulation m the Pioneer 
11 data is just as large at 2.42 Rs as in the 2.65 Rs gap, as determined from the proton 
profiles published by Van Allen et al. and by Fillius et al. [1980], where data from the latter 

can be seen in our Figure 2. 

Thus the available anisotropy data do not address the relative merits of aimpoints m the 
three proton gaps, but the Voyager data in Figure 11 do indicate more definitively that 3.2 
Rs is in a good region from the standpoint of high anisotropy. The higher proton intensity at 
2 65 Rs suggests that one might find a larger anisotropy there, but the published Pioneer 11 
show no discernable difference from that found elsewhere, including at 2.8 R s where Figure 
11 shows a low anisotropy level. Additional input from the Pioneer 11 investigators (e.g., 
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Van Allen et aL, Fiffius ti al.) would be desirable to allow study of the available, unpublted 
^Thegenerd 3^ Rs is conceivably related to dust abocrption 

is the more likely explanation. This needs more work and might also reveal something 
interesting about Urge particle distributions in the main rings which emit the parent neutrons 
for CRAM) Until one understands this better, I would be inclined to move the most critical 
^ g SO^aimpoL as far away from Enceladus’s orbit as possible and would again choose 
2 7 o’ver 3:2 R S The available anisotropy data do not rule out an aimpoint even at 2.42 R s . 

FOt Are a t^ a !ub°MeV ion and electron measurements inside 5 R s relevant for constraints 
on E-rfng densities? Considering the not-so-pathological scenario that I have conc«-d for 
the Vovaeer 2 LECP response to MeV electrons in this region, I will personally ha 
confMence^at all in these data until laboratory experiment, (easily done) estabbsh that 

internal back, mattering of MeV electrons cannot explain the obsM 7‘““' 

that Van Allen’s early reservations about instrument backgrounds from the electrons, having 

a'drift resonance peak in differential spectra near 10-1.2 MeV (i.e„ not sufficren o 

the LECP sun shield used to check for background response), require much further scru y 

wuh faboratory and theoretical investigations before final decisions are made about E-rmg 

h “ The microsignature observations may have similar problems ’ “ V ^rECr^p»thS 
satellite orbits. 
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